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Abstract. We develop a model of production of the very dense clumps of DM in RD epoch 
due to the accretion of DM on the loops of cosmic strings as the seeds. At some time the 
loops disappear, for example due to the gravitational radiation, and the remaining dense 
clumps produce the enhancement of the annihilation signal. We take into account the ve- 
i-G i locity distribution of the strings, and consider the two extreme regimes of DM annihilation: 

*i i fast decay and continuous evaporation. The produced annihilation flux of gamma radia- 

tion is detectable, and for some parameters of DM particles and the strings can exceed the 

+r* ■ extragalactic flux of the gamma-radiation observed by Fermi. For the fixed parameters of 

DM particles (e.g. neutralino with fixed masses and cross-section of annihilation) one can 



obtain the limits on the basic string parameter, tension if, which is stronger than (more 
CN ■ general) limits obtained from WMAP observations, cosmological nucleosynthesis and gravi- 

tational lensing. In particular for the neutralino with 100 GeV mass we exclude the interval 
IT) I 5 x 10" 10 < Gfi/c 2 < 5.1 x 10" 9 . 
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1 Introduction 

The linear topological defects — cosmic strings can be formed in the early cosmological phase 
transitions (see for a review [1], [2]). Along with the infinite strings there is possibility of 
closed loops formation in the network of curved cosmic strings due to their interconnections. 
According to numerical simulations, after a long transient stage a true scaling regime sus- 
tained, then the typical distances between the strings and the coherence length both scale in 
proportion with the horizon scale [3]. A string loop formed at the cosmological time ti has 
the length I ~ acti, where in the scaling regime a ~ 0.1 according to [3] and [4], although 
other values of a were obtained in other works (see for example [5], [6]), down toa~ 1CP 3 . 

A fundamental characteristic of the string is the mass per unit length \x = Mi J I or the 
tension, which is of the order of symmetry breaking energy squared rj 2 . For example, the 
grand-unification-scale strings have Gp/c 2 ~ 10~ 6 , where G is the Newton's constant. There 
are several restrictions on fi. From CMB observations it follows G/i/c 2 < 2 x 10~ 7 [7]. The 
bound G/x/c 2 < 10~ 7 was obtained from the study of nucleosynthesis [3]. Search for the pairs 
of galaxies images consistent with the gravitational lensing of the cosmic string presents the 
limit Gfj,/c 2 < 3 x 10 -7 at 95% confidence level [8] . In [9] the formation of stars in the first dark 
matter DM haloes seeded by the loops was considered. It was found that Gp/c 2 < 3 x 10~ 8 to 
avoid collision with WMAP data on the reionization redshift. Searches for the gravitational 
wave bursts from strings by LIGO provide the joint constraints on the strings parameters 
(/i and interconnection probability) [10], but Gp/c 2 is only weakly restricted in comparison 
with the constraints above. And finally, the strongest bound Gfi/c 2 < 4 x 10 -9 was obtained 
from the pulsar timing [11]. 

In this work we present new constraint on \x which was obtained from the DM particles 
annihilation in the dense clumps seeded by the loops at the cosmological stage of radiation 
dominance. The direct detections of DM particles are the promising but still elusive ex- 
perimental problems, therefore the search of indirect signature of the DM is important for 
clarifying the DM origin. The promising indirect signature -- DM particles annihilation 
would proceed more efficiently (it can be boosted by several orders) if the Galactic halo is 
filled by the dense DM substructures or DM clumps. The cosmic string loops produce very 
dense clumps due to their early formation. Only low velocity loops can produce the clumps. 
The probability of the low velocity loop formation is very small, but even tiny fraction of the 
formed loops may produce the dense clump population and significant annihilation signal. 

The clumps formation at the radiation dominated (RD) stage was studied in details by 
[12]. In the particular case of loops' density perturbations the maximum density of clump is 
restricted due to adiabatic expansion of the already formed clump after the loop gravitational 
evaporation. We found the modification of this restriction in the case then the loop decays 
before the clump virialization. In this case the clumps can reach density p c \ ^> 140p cq , 
where p cq is the density at equality. The comparison of the resulting annihilation signal 
with the Fermi-LAT data allows us to obtain the restriction on the string parameter fi. 
It must be pointed out that our constraints were obtained by supposing that the DM can 
annihilate, and we take the ~ 100 GeV neutralino as the most promising particle candidate. 
The constraints will be different for other DM models. In this sense the obtained constraints 
must be considered as joint constraints on the properties of the string loops and DM particles. 



2 Initial speed of the loops and rocket effect 

Here we consider the influence of the initial velocities of the loops and rocket effect on the 
evolution of perturbations and clumps formation. 

The necessary condition for the clump formation is the low velocity of the loop [12]. 
Only results for average initial velocity of formed loops were presented in literature but the 
distribution over velocities is possible. We interested in low velocity end of the distribution, 
because the clump forms only if the seed loop stays near the center of the clump during its 
evolution. Loops can be formed by intersecting long string segments or by self-intersection of 
long strings. We suppose that the probability of velocity components of the loops is simply 
Gaussian with mean value at the correlation length scale (vf) ' — 0.15c [13], and therefore 
the probability of full initial velocity is 

Even if the process of a loop's formation involves the intersection of two strings with large 

velocities, it does not necessarily means that the resulting velocity will be high. 

The displacement of the loop beginning from its birth moment ij till the full decay 

moment td is 

t& 

v(t)dt 



Ar = a ^Jiw' (2 ' 2) 

u 

where the peculiar velocity is v(t) = Via(tj) / a(t) . We require that the displacement Ar is 
smaller in comparison with the loop's radius l/(2ir). From this condition and (2.2) we obtain 
the restriction on the loop's velocity 

v i t i ln(t d /t i )<l/(27r). (2.3) 

For the probable parameters of the strings td/U ~2x 10 5 and the dependence in (2.3) is only 
logarithmic. Therefore, by using (2.1) we can estimate the probability of the low velocity 
loop formation as 

V^ = .^. (-) 

As we will show below even the tiny fraction (2.4) of the formed loops may produce superdense 
clumps and observable annihilation signals. 

Now we consider the rocket effect. Velocity of the loop grows linearly with time v r = 
3TpG/j,t/(5l), where Tp ~ 10 [1]. Turnaround moment of the clump corresponds to t\ A ~ 
500£ 2 , and the relative displacement of the loop during clump formation is 

- I v r dt ~ 1.5 x 10 - V_ 8 o^ < 1, (2.5) 

U 

where ^_g = G/i/(10 _8 c 2 ). Therefore, for the small loops formed at radiation era the large 
rocket displacements are not achieved. 



3 Evolution of clumps around evaporating loops 

We solve the same equation as eq. (2.7) in [12] 

_,d 2 b |\ 3 ] db 1 Tl + * ,' 
X{x + 1) d^ + l+ 2 X Tx + 2 —~ h 



0, (3.1) 



where x = a{t)/a eCL is used as independent variable, r = ab{x)^ is the physical radius of the 
spherical shell and £ is its comoving coordinate, $ is the density perturbation of DM Sp/p 
inside the spherical shell. This equation describes the evolution of the clump's radius r in 
terms of the function b(x). The only quantity one need to modify is the <3?. In difference with 
[12] we allow the dependence of 3> on the time: steady decrease in the continuous evaporation 
approximation and step-like in the fast decay approximation. The evolution of clump stops 
when dr/dt = or equivalently db/dx = —b/x [12]. The density and the radius of the clump 
at the moment of its maximum expansion are 

3M N 1/3 



Pmax — Peq^max^max' "max — I . I ; \"-^J 

V 4vrp max / 

where 6 max and x max are the values at the moment of the stop. After the turnaround the 
clump virializes by contracting twice in radius, and the resulting density increases by factor 
8 in comparison with (3.2). 

Let us name the spherical region with a volume (4-7r/3) • (//2-/r) 3 as a "string volume". 
Then the fraction of string mass Mi = pi to the mass Mf^ M of DM inside the string volume 
at the moment of the string birth ti = l/(ac) is simply 
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M p 



(3.3) 

t=U 



where Mr = 1.6 x 10 3 /^LgCKg \M & . The (3.3) was calculated from the law of DM density 
evolution ex a -3 according to the known solutions of the Friedmann equations. The fraction 
(3.3) gives also the value $ of the density perturbation inside the string volume at the 
moment of string birth ij, and in the most interesting cases <3? S> 1, with the production 
of the superdense clumps. The strings with Mi = Mf) M (t = ti) born at the time tp = 
3.9 x 10 _6 /x?_ 8 a|^ 1 teq {xp = 2 x 10 -3 /i_8«o\)- We consider only the most dense central part 
of the clump inside the string volume, where the annihilation proceeds most effectively. This 
central region of the clump can be refereed as a clump core. The oscillation of string itself 
does not distort the core because the linear segment of the string has small gravitational 
potential. The only common (mean over the string volume) potential of the loop attracts 
the DM. The outer regions of the clump form through the secondary accretion of DM and 
have the density profile p(r) oc r~ 9 ' 4 at the sufficiently high distance from the center of the 
clump. Therefore the annihilation concentrates near the clump core. 

4 Continuous evaporation and fast decay approximations 

The characteristic loop lifetime due to gravitational waves emission is r ~ lc/(GpT), where 
r ~ 50 is a numerical coefficient [14] . One possibility is that the loop losses mass continuously 
according to mean equation dMi/dt = —TGp 2 /c, but the more reliable approach is to assume 



the sudden decay after the time interval r from the birth moment t{. At t > ti + r the loop's 
configuration substantially changed, so the loop is likely to self intersect. The resulting 
daughter loops will fly away at high speeds. We use the last approximation henceforth 
as the main but also present the results for the continuous evaporation approximation for 
comparison. 

DM clumps formation at the RD stage was explored in [12]. In the particular case for the 
clumps which are seeded by loops of cosmic strings one have $ ~ Mi/M, but in many cases 
the maximum density is only p c \ ~ 140/> e q due to adiabatic expansion of already formed 
clump during the seed loop gravitational evaporation [12]. Really, the universal Poincare 
adiabatic invariant conservation J = §Y1 Pidqi = const for the clump with additional loop 
mass inside implies M tot i? = const or p c \ oc M t ~ t , where M tot = Mi(t) + Mdm is the total 
mass of the loop and DM. For the constant mass loop the clump forms with the density 
p c i ~ 140 / o eq (M//MDM) 3 (see eq. (3.4) in [12])). After the subsequent loop decay the density 
lowered due to adiabatic invariant in proportion ~ (Mj^m/Mi) 3 till the value p c \ ~ 140/3 cq . 

We argue that the above mentioned argument of the adiabatic invariant conservation 
is not applicable, if the loop decay occurs before turnaround moment (detachment from the 
cosmological expansion and the clump virialization) . Really, in this case DM particles move 
not at the orbits around the loop but along the radial trajectories. The loop's decay leads 
only to the change of the particles acceleration. The evolution of clump slows down but 
continues under the influence of velocities db/dt obtained before the decay. We can estimate 
the processes by the following manner. For the clumps under consideration the conditions 
i< 1, $>1 and <£>£ <C 1 are valid almost all the time before turnaround. Initially the 
evolution goes due to large value of $ and the initial velocities db/dx at t = U are not 
important. At this stage one can neglect the first term in (3.1) and the approximate solution 
is b ~ 1 — x<3>/2 [12]. If the turnaround occurs before tj the moment of turnaround can be 
estimated as xta ~ 1/$ [12]. In the opposite case Xd < xta just after loop decay (at x = x<j) 
we must put $ = in (3.1) and the velocity at this moment db/dx = — <J>/2 becomes greater 
then the last term in (3.1). At x > Xd we can neglect the last term but leave the first one. 
This leads to the red-shifting of velocity as 



db _ ®x d 
dx 2x 



(4.1) 



and to the corresponding logarithmic decrease of b. From the condition db/dx = —b/x we 
obtain the new turnaround moment: 

XT ^ Xd ^{^A). (4.2) 

We see that at the sufficiently small $x^ <C 1 the clump may not forms at all, because the 
xta will be exponentially large. For the moderately small values §Xd <C 1 the clump forms 
but with small density. The value &Xd can be expressed as 

<bx d a 0.9^y 8 2 a-? /2 r 5 - 1/2 . (4.3) 

For the constant mass loop xta ~ 1/^ [12] and (4.3) is ~ x^/xta- This value is close to 
unity at /x_g ~ 1, so we expect the change in the character of the clump formation process 
near /i_g ~ 1. 



In the continuous evaporation approximation the rate of loop mass evaporation due to 
gravitational radiation is dMi/dt = —TGp 2 /c. After integration we have 



«o.i 



(4.4) 



where the M/(t,) = pati is the initial string mass at the moment of its birth tj. One need 
to generalize the evolving seed mass as the fuse of fluctuation growth. The only quantity 
one need to change is $. From the known solutions of the Fridmann equations one has the 
dependence t = t{x) and by using (3.3) we finally obtain 

H X , XI ) = 2XW ^^ M ^ X 1. (4.5) 

Xi Mi{t{xi)) 

This expression is valid for <I> > 0. If formally $ < we put <I> = in (3.1). This means 
that the string had totally evaporated (its mass is zero) and the subsequent clump evolution 
proceeds only under DM self gravity and due to inward velocity boost which appeared before 
the string decay. The velocity boost leads to the perturbation growth even after full evapo- 
ration of the seed loop. The clump virializes with some density /? c i(£ta)- If the string mass 
goes to zero after the turnaround the adiabatic expansion of the clump occurs only due to 
the loop mass remnant M^(£ta)) and the resulting clump density is 

u v ( M DM \ Pci(*ta) Ia r s 

Pcl = Pd(iTA V^TA) + M DM J = (1 + *(S TA ))3- (4 ' 6) 

In the case td > £ta this density is greater in comparison with density p c \ ~ 140p cq in the 
fast decay approximation. 

5 Numerical results 

We solve Eq. (3.1) numerically in the two above approximations. In the approximation of fast 
decay the loop decay at the moment td = U + r. This means that at t < td the mass of the 
string is constant, but at t > td the string had totally disappeared (its mass is zero), and the 
subsequent clump evolution proceeds only under DM self gravity and due to inward velocity 
boost which appeared before the string decay. The velocity boost leads to the perturbation 
growth even after full evaporation of the seed loop. 

We consider only the most dense central region of the clumps, which gives the main 
contribution to the annihilation signal. These are regions inside the string volumes. As the 
first approximation we consider these regions as homogeneous. We find the density of the 
clump in dependence of Xi and p. The turnaround moment is calculated numerically from 
the condition db/dx = —b/x and the solution of (3.1). Clumps density is obtained according 
to (3.2) and (4.6). If the turnaround moment precedes the loop decay, we put the resulting 
clump density p = 140p cq according to adiabatic invariant argument conservation of [12]. 

The results of calculations for clumps density in the fast decay approximation are shown 
at the figure 1. As it was expected from (4.3), the condition xta — %d is satisfied near 
fi-8 ~ 1) and the regime of clump formation changes near p^$ ~ 1 because at larger /i_g 
the turnaround occurs before the loop decay. The similar figure can be presented for the 
continuous evaporation approximation, but with smoother surface break and with greater 
density of the clumps. 
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Figure 1. Clump density p in the units of density at matter-radiation equality p c(i in dependence 
on the loop birth moment Xi and parameter /i_g = G/i/(10~ 8 c 2 ). The break of the surface down to 
value p = 140p C q corresponds to the proximity of turnaround and loop decay moments. 



6 Loops and clumps distributions 

The length distribution of cosmic strings' loops in the interconnecting network was obtained 
in [9] in the form 

Ndl 

dnio °P ~~ c 3/2 f 3/2/5/2 ' t 6 - 1 ) 

where N ~ 2. The evaporating mass cutoff must be superimposed on the distribution (6.1) 
at the every particular time. If we neglect (temporary) the loop evaporation, then the mass 
fraction of the universe in the form of loops at the time t eq is 



^M = . 42^_/ 8 2 

Pcq 



M 



-3/2 



dMi 



(6.2) 



In terms of cosmological density of clumps (a fraction of DM mass in the form of clumps) 
the distribution (6.2) (by using (3.3)) can be rewritten as: 
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(6.3) 



where P\ v is given by (2.4). Strings decay but the clumps survive, therefore there is no need 
to cut of the clumps mass spectrum at the string evaporation scale, and the (6.3) is the 
real distribution of clumps at MD epoch. The density of these clumps was calculated in the 
Section 5. 
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Figure 2. Upper limits on (crv) (in units 10 -26 cm 3 s" 1 ) in dependence of the string parameter 
/i_8 = Gfi/(10~ s c 2 ). The solid lines show the limits for the masses of DM particles (from up to 
down) m x = 10 TcV, 1 TcV, 100 GcV and 10 GcV in the fast decay approximation. The limits were 
obtained from the comparison of the calculated signals and the Fcrmi-LAT data. The upper and lower 
horizontal dashed lines show the typical and minimal possible cross-section values, respectively. The 
dotted lines show the upper limits in the continuous evaporation approximation for the same masses. 



The low mass cut of of the clumps distribution is determined by the process of kinetic 
decoupling of the DM particles. At earlier times the DM particles strongly frozen in the 
radiation and do not move toward the loop. In contrast to the ordinary inflationary density 
perturbations the diffusion and free streaming effects are not important for the minimum 
mass of the clumps. This is because the forming clump subjected mainly by the strong 
gravitational pull of the central loop and evolve nonlinearly long before the equality moment 
t e q- The kinetic decoupling temperature for ordinary neutralino weakly depends on the 
particle mass T^ oc m x and, for example, for toioo = m x /(100 GeV) = 1 and for typical 
SUSY parameters T^ ~ 25 GeV with corresponding cosmological time td ~ 1.2 x 1CU 3 s. The 
loops which formed at the moment td have masses Mi 



r 3/2 



— 1 /2 

2.5 x 10 - 'm 10 Q /i_gao.i-^0 and 



K™JM t 



1/2 



100 

2 x lO- 15 m- 3 o /4 ^ lM0 



the minimum clump's mass is therefore M c i jm ; n 

according to (3.3). These minimum mass clumps can reach densities p c \ ~ 3 x 10~ 4 g cm" 3 

if /i_g ~ 0.4 (see figure 1). 



7 DM annihilation 



The clumps under consideration have very large densities and the gamma-ray flux from DM 
annihilation inside the clumps may exceed the observational limits for some values of string 
parameter /i_8- Let us consider the neutralino (the most popular DM candidate) annihilation 



R 

in the clumps. Annihilation rate of neutralino in a single clump -/V ann = 2r) n o4ir(av) J n x r 2 dr, 

o 
where 77^0 ~ 10 is the neutral pion multiplicity, n x is the number density of particles inside 

clump, R ~ (3M/47T/)) 1 ' 3 and (av) is the annihilation cross-section (averaged product with 

velocity). We consider the annihilation channel with 7r° productions and decays 7r° — > 27. 

The cumulative gamma-ray signal from the clumps in the angular direction ip with respect 

to Galactic center can be expressed as 

M E > m ,„/2,« = 1.9 x 10-- (^r 10 -J'^,.- 1 W)U. (7.1) 

where 

W))ao = [d& ( ** ,) [ -4- ( PH(r) ,) , (7.2) 

\ vvj/mi j ^ cl Vo.3 GeVcm" 3 / 7 8.5 kpc Vo.3 GeV cm" 3 ; ' V ; 

l.o.s. 

and the last integration goes along the line of sight. For the halo density profile pn(r) we 
use the NFW profile [15] with the scale a = 20 kpc, halo mass M^ = 10 12 M Q and virial 
radius Rh = 200 kpc. The lower mass limits in the integration d£ c \ was estimated in the 
previous section. This limit weakly depends on m x through the Td{m x ) dependence. We 
took M/ jinax ~ 1.6 x lO 3 /i 3 _ 8 M0 (this seed mass corresponds to the clump's formation time 
near t eq ) as the upper limits of the integration, and the dependence of the final result on 
-M"z,max is weak. 

We compare the calculated signals with Fermi-LAT diffuse extragalactic gamma-ray 
background J \, S (E > m n o/2) = 1.8 x 10~ 5 cm" 2 s" 1 sr" 1 [16]. To obtain the most conser- 
vative limit we compare J b s with the calculated signal in the anti-center direction ip = ir. It 
gives the upper limit on (av) in dependence of fi-s = G/j,/(W~ 8 c 2 ). The results are shown 
at figure 2. We consider the several values of the neutralino mass: m x = 10 TeV, 1 TeV, 
100 GeV and 10 GeV. The mass m x influences the result mainly through the factor m~ 2 un- 
der the integral (7.1) and through the low mass limit in the loops distribution, which weakly 
depends om m x . 

For example, for the mass m x = 100 GeV in the case of typical neutralino cross-section 
(av) ~ 3 x 10 -26 cm 3 s _1 (this value corresponds to the thermal production of DM particles) 
the limit excludes the range of parameters 0.05 < /i_g < 0.51 in fast decay approximation 
and 0.1 < ;U_8 < 1.16 in the continuous evaporation approximation. If we take the minimal 
allowed value (av) = 1.7 x 10~ 30 m~" 00 cm 3 s~ 1 [17] and m x = 100 GeV the excluded regions 
are 0.16 < /i_s < 0.43 and 0.27 < /i_s < 1.07 in the same two approximations. 

8 Conclusions 

In this work we explore the possibility of dark matter annihilation in the dense cores of 
clumps that were seeded by cosmic string loops before the matter-radiation equality. We 
calculate the evolution of the clumps around evaporating loops. Only the low-velocity loops 
result in the clumps formation. At the same time, even the low-velocity tail of the loop 
distribution produces DM clumps with the observable signature in the annihilation products. 
The adiabatic argument conservation doesn't prevent the formation of clumps with densities 
p c i S> 140/Jeq, if the decay of the loops occurs before the time of the clumps virialization. 
Therefore the DM clumps produced by the loops can be the very dense objects with high 
cumulative luminosity in the gamma-rays. 



The combined constraints on the loops parameters (/i and the distribution over lengths) 
and parameters of DM particles were obtained. For the 100 GeV neutralino DM the range 
5 x 10~ 10 < Gfi/c 2 < 5.1 x 10~ 9 was excluded because of the huge gamma-ray annihilation 
signal above the Fermi-LAT data for these parameters. Along with the preferable neutralino 
dark matter candidate with the mass 100 GeV, we consider the masses 10 TeV, 1 TeV, and 
10 GeV and explore the two limiting assumptions about the loop evolution: fast decay and 
continuous evaporation approximations. The corresponding constraints on the annihilational 
cross-section are shown at the figure 2. We restrict ourself only by the preferable value 
a = 0.1. The results will change for different a. 

In the case of the cosmic superstrings, the reconnection probability can be <C 1, resulting 
in the much higher number density of loops. This could lead to even stronger constraints in 
comparison with the presented in this paper. 

The Fermi-LAT data are used as the upper limit. In principle the annihilation of DM 
in the clumps can explain the observed signal for the particular values, for example G^i/c 2 ~ 
5 x 10 -10 . The necessity of the DM annihilation can arise if the ordinary astrophysical sources 
give too small signal in comparison with the observations. 
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